Brownian motion exhibiting absolute negative mobility 



Ralf Eichhorn, Peter Reimann, and Peter Hanggi 
Universitat Augsburg, Theoretische Physik I, Universitdtsstr. 1, D-86135 Augsburg, Germany 



o : 
o . 

(N ■ 

?h ; 
<: 

t> ■ 



43 

o 



CO 

o3 



i 

C 

O 
o 



> 

m 
o 

O 



-a 
c 

o 
o 



X 



We consider a single Brownian particle in a spatially sym- 
metric, periodic system far from thermal equilibrium. This 
setup can be readily realized experimentally. Upon applica- 
tion of an external static force F, the average particle velocity 
is negative for F > and positive for F < (absolute negative 
mobility) . 

PACS: 05.40.-a, 05.60-k, 02.50Ey 

When a system at rest is perturbed by a static force, 
we expect that it responds by moving into the direction 
of that force. The rather surprising opposite behavior 
in the form of a permanent motion against a (not too 
large) static force of whatever direction is called abso- 
lute negative mobility (ANM) If the unperturbed 
system is at thermal equilibrium then ANM is impossi- 
ble since it could be exploited to construct a perpetuum 
mobile of the second kind. Legendary but quite com- 
plex non-equilibrium systems which do exhibit ANM are 
"donkeys" B . In this Letter we demonstrate that even a 
simple, structureless Brownian particle can exhibit ANM 
away from thermal equilibrium in a setup which can be 
easily realized experimentally — a scenario hitherto com- 
monly considered as impossible. 

Going in medias res, let us consider the overdamped 
2-dim. Brownian motion 



iji(t) = -d x V(x{t),y(t)) 
r)y(t) = -d y V(x(t),y(t)) 



F 



(1) 



where r\ is the viscous friction coefficient and V(x,y) 
is the hard-wall potential from Fig. la, confining the 
Brownian motion to a "corridor" along the y-axis with 
symmetric, periodic obstacles. Further, £ x (t) are thermal 
fluctuations, modeled by unbiased Gaussian white noise 
with (£x{t) £x(s)) — 2i] ksT S(t — s), where ks denotes 
Boltzmann's constant and T the temperature. If £ y (t) is 
an independent, second thermal white noise source then 
(Q) is an equilibrium system and the average particle cur- 
rent 



(y) ■= 



lim 

t — >oo 



y(t) - y{t ) 
t-t 



(2) 



always runs into the direction of the static force F. Fig. 
lb depicts the corresponding (y)--F-characteristics for the 
simplest non equilibrium model, namely a symmetric di- 
chotomous noise £ y (t), with the promised ANM as its 
most outstanding feature. 

Our first remark is that the so-called ratchet effect ||] 
is characterized by a current (y) which is non-zero for 
F = and does not change its direction within an entire 
neighborhood of F — 0. This effect thus inevitably in- 
volves some kind of asymmetry (for F = 0), whereas our 



present system is perfectly symmetric. Second, so-called 
differential negative mobility (or resistance) pj is typified 
by a negative slope of the (y)-i ? -characteristics away from 
F = 0. Thus, in both cases the salient feature of ANM 
is absent, namely a current (y) which is always opposite 
to the (not too large) force F, independently of whether 
F is positive or negative. Finally, we mention that ANM 
has also been observed in semiconductor devices [jlj and 
in models for coupled Brownian motors However, 
in the first case it has an entirely quantum mechanical 
origin and in the second case it is a genuine collective 
effect, without leaving room for any kind of classical, 
single-particle counterpart. Accordingly, the respective 
physical mechanisms are completely different from ours. 

Returning to the ANM in Fig. lb, its origin can be un- 
derstood as follows: Consider a (moderately large) time 
interval t during which the dichotomous noise £ v (t) in 
(|]) is constant and F tot :— C y (t) + F > 0. When starting 
in one of the "corners" between the right "corridor wall" 
and any of the adjacent obstacles (see Fig. la), the par- 
ticle first closely follows the right "corridor wall" until 
it hits the next obstacle, then "slides down the back" of 
that obstacle, and afterwards performs a "free fall" in the 
y-direction. Since the lateral extension of the obstacles b 
exceeds half the corridor width B/2 the particle then hits 
with a high probability q the next obstacle on its way and 
ends up by being trapped in the corresponding (left) "cor- 
ner" . In order to avoid this trap, the particle has to ther- 
mally diffuse at least over a distance b—(B — b) = 2b — B 
in the positive ir-direction during its "free fall" in the y- 
direction. With increasing force F to t, the available time 
and therefore the probability p := 1 — q of such a diffusive 
displacement decreases, implying that the particle travels 
on the average a shorter distance along the y-axis during 
the time r. Since an analogous consideration applies for 
time-intervals r with F tot = £ y (t) + F < 0, the particle 
motion on the average acquires a bias in the direction 
opposite to the static force F, i.e. it exhibits ANM. 

In order to quantify this argument, we first note that 
the above mentioned probability p of avoiding a trap (for 
Ftot > 0) can be approximated as B 



, 1 1 / 2b -B 

p{Ftot) ^2-r ii {v2 = mT 



(3) 



where erf(x) := 27r~ 1 / 2 J* e~ y dy. With probability p, 
a particle thus covers in addition to the "basic distance" 
of approximately 3L/2 another period L (see Fig. la). 
It then avoids the second trap on its way with approx- 
imately the same (relative) probability p as in (H), i.e. 
a second period L is covered with (absolute) probabil- 
ity p 2 etc., see Fig. la. If the maximal traveling distance 



1 



(avoiding all traps) is of the form (3/2+N)L with N e N, 
the average traveling distance Ay(r, F to t) thus follows as 
L[3/2+p + p 2 + ...+p N ]. Neglecting that the "free trav- 
eling speed" v y :— Ftot/v is slightly reduced when the 
particle "slides down the back" of an obstacle, we obtain 
(3/2 + N)L = v v t and hence 



Ay(T,F to t)=L 




[p(Fu 



l-p(Ftot) 



(4) 



This expression remains a decent interpolation even if v y r 
is not of the form (3/2+N)L. Symmetrically, for F to t < 
the average traveling distance is — Ay(r, —F to t), implying 
for the net average current (0) the approximation 



(y) 



r °° dr p(r) [Ay (r, A + F) - Ay (r, A - F)} 
2j oo drp(T)T 



(5) 



where ±A are the two states of the dichotomous noise 



£ y (t) and p(r) 



7e" 



-7T 



is the distribution of sojourn 



times, i.e. 7 is the flip rate. 

The agreement of our analytic prediction (|^) with the 
simulations in Fig. lb is remarkably good in view of the 
various underlying approximations. In particular, our 
assumption that the particle covers at least a distance 
3.L/2 during the time r renders doubtful unless v v t — 
FtotT/rj > 3L/2. To fulfill this condition for all the forces 
Ftot — A ± F notably contributing in ([5]) requires that 
A - \F\ > 3ji]L/2. Indeed, for A < 3jr)L/2 ANM is 
found to disappear in numerical simulations of (|l|). 

On the other hand, ANM is expected to subsist for 
numerous generalizations of our original model ([!]). Fig. 
2 exemplifies a setup with a 2-dim. array of obstacles. 
For symmetry reasons, the current (||) remains exactly 
the same as in Fig. 1, but the parallelization now ad- 
mits to simultaneously transport a much larger number 
of particles. Such a device can be readily realized by a 
modification of those studied experimentally in |tJ and 
theoretically in J§| . We emphasize that while these mod- 
ifications of the experimental setups are straightforward, 
the physics, however, is completely different. A further 
step towards a realistic experimental system is achieved 
by choosing 



ty(t)=tith(t) + f(t) 



(6) 



where t;th(t) is another thermal white noise like £ x (t) (but 
statistically independent) and f(t) switches periodically 
between ±A with p(r) = S(t — r ac ). For not too weak 
driving f(t) and with the appropriately adjusted defini- 
tion F to t '■= f(t) + F, the corrections in (0), (jij) due to 
the thermal noise £th (t) are small and thus (|5j) remains a 
valid approximation provided A — \F\ > 3r/L/2r ac . 

An experimental realization of the system (Q), @ 
along the lines of Q,^] is presently under construction in 
the labs of C. Bechinger and P. Leiderer. Henceforth, we 
focus on such experimentally realistic parameter values 
in our quantitative examples, see Fig. 3. In particular, 



the agreement of (||) with the numerical simulations in 
Fig. 3a for r QC = 1 s is again rather good in the parameter 
range LF I < 0. 14 pN compatible with A- \F\ > 3?/L/2r ac . 

In (Jfj), we have completely neglected the possibility 
that a trapped particle may escape from the trap due to 
the ambient thermal noise. This is justified as long as r 
is much smaller than the mean escape time T esc (F to t) out 
of a trap. Turning to the opposite case r ^> T esc (F to t), 
we start by calculating the time which a particle needs 
to advance by one period L along the y-axis for F to t > 0: 
This time is approximately L/v y if the trap within such 
a period is avoided and L/v y + T esc (F to t) otherwise. The 
respective probabilities are p and 1 — p, approximated 
by (||), i.e. the average time to cover one period L is 
pL/v y + [l- p] [L/vy + T es c(Ftot)] ■ The resulting average 
traveling distance during the (large) time r 3> T esc (F to t) 



Ay(T,F tot ) 



tL 



^- t +[l~p(Ftot)}T esc (F tot ) 



(7) 



Thus, if those large r dominate, Ay /r becomes indepen- 
dent of t and ([5]) can be rewritten as 



(y) = - 2 



Ay(r,A + F) Ay(r,A-F) 



(8) 



independent of p(r). For small F to t > 0, the first term in 
the denominator of (0) dominates and hence Ay increases 
in the expected linear response manner with increasing 
Ftot- As Ftot gets larger, 1 — p approaches 1 (cf. (||)) and 
the escape time r esc increases very fast (cf. (||) below), 
implying the existence of a maximum and a subsequent 
decay of Ay. As a consequence of this increasing "stick- 
iness" of the traps with increasing F tot , we recover 
once again ANM in (^) provided A is sufficiently large. 

Focusing on the model (|l|), (||), one can approximate 
Tesc(Ftot) by the mean first passage time from x = to 
x = b/ sm.9 of the auxiliary 1-dim. dynamics rjx(t) = 
—Ftot cos 9 + £, x (t) with a reflecting boundary at x = 0, 
reading 110] 



b 2 77 e a — a — 1 



ksT a 2 sin 2 1 



a : = b F^t cot 9 /k B T ■ 



(9) 



The agreement of with the numerical simulations 

in Fig. 3a for r ac = 25 s is quite satisfactory. In particu- 
lar, the predicted ANM is indeed recovered. 

We emphasize that the basic physical origins of ANM 
are completely different in the small- and large-r regimes 
as quantified by (^) , (||) and (0) , (||) , respectively: In the 
former case, escapes out of the traps are negligible, while 
transient "first-trapping events" after each jump of Ftot 
provide the crucial mechanism for ANM. In the latter 
case, these transients are negligible, while the "re-escape 
events" are now at the origin of ANM. This remark- 
able feature that two completely different physical mech- 



2 



anisms both support one and the same phenomenon, sug- 
gests that ANM will also be present in the so far disre- 
garded intermediate-r regime. Furthermore, on the ba- 
sis of our above physical insight an immediate educated 
guess is to add up (|j) and (0) and then evaluate (||). 
Both these predictions are nicely confirmed by Fig. 3b. 

A more sophisticated analysis can again be based on 
our usual assumption that a particle always closely passes 
by the leftmost edge of any obstacle attached to the 
"right corridor wall" when F to t > 0, as indicated in Fig. 
la. Consequently, the traveling times through any period 
L are governed by one and the same probability distri- 
bution tp(t), independent of the particle's past (Markov 
property). Similarly as in (Q), this distribution is approx- 
imately given by 



i>(t)= P 5(t-n) + (i-p)e(t-T 1 y 



-({-Tl)/T e 



(10) 



where &(t) is the Heaviside function and t% := L/v y . In 
this way, the original 2-dim. problem can be approxi- 
mately reduced to a 1-dim., uni-directional hopping pro- 
cess characterized by if>(t). Such processes have been 
analyzed in detail in the context of renewal theory [|ll| . 
Along these lines, we obtain for the Laplace transformed 
displacement Ay(s, F tot ) := / °° dt Ay(t, F tot ) e~ ts there- 
suit 



Ay(s,F tot ) = ±-^ 
s 1 



1 



xl>(a) 



TlS 



where ip{s) is the Laplace transform of ip(t). While 
the first two factors on the right hand side of ([ll]) are 
well known jll], the last factor accounts for the fact 
that the particle actually proceeds continuously rather 
than in discrete jumps of length L. Moreover, after the 
Laplace back-transformation of (^), a final transforma- 
tion Ay{T,F tot ) h-> Ay(r - 3L/2v y ,F tot ) + 3L/2 is re- 
quired since the "basic distance" 3L/2, which the parti- 
cle covers before encountering the first trap (see Fig. la 
and below (0)), is not yet taken into account by ([ll]). For 
very small and large r- values one then recovers our previ- 
ous results (^) and ([?]), respectively, while for more gen- 
eral r- values, a numerical evaluation of the Laplace back- 
transformation is necessary. A more detailed derivation 
of these analytical results will be presented elsewhere. A 
typical example is depicted in Fig. 3b, in good agreement 
with the numerical simulations. 

In conclusion, we have demonstrated that a single, clas- 
sical Brownian particle in a periodic, symmetric 2-dim. 
potential landscape can exhibit the prima facie quite as- 
tonishing phenomenon of absolute negative mobility un- 
der suitable far from equilibrium conditions. In general, 
the effect is simultaneously supported by two completely 
different physical mechanisms and, in contrast to jjj, is 
not restricted to adiabatically slow non-equilibrium per- 
turbations. The phenomenon is moreover robust against 



modifications of the potential. It can occur in practically 
any potential landscape that provides "traps" with in- 
creasing "stickiness" as external force strengths increase. 
The setup discussed here is particularly suitable for an 
experimental realization along the lines of Jt|,||J^] ■ 
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HA1517/13-4 and the Graduiertenkolleg GRK283. 
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FIG. 1. (a): Hard-wall potential V(x,y), denned as zero in 
the white and infinity in the black regions. It confines the mo- 
tion of a particle to a "corridor" of width B along the {/-axis 
(the white regions outside this corridor are of no interest). 
The black obstacles are arranged in a symmetric, periodic 
manner (180° rotation symmetry) and are characterized by 
the parameters 9 and b (their extension in {/-direction is ne- 
glected). Note that b > B/2. Several typical traveling routes 
of a particle under the action of a static force Ftot > along 
the y-axis are indicated as arrows together with their prob- 
abilities (q := 1 — p). (b): The current (^) versus the static 
force F in (|l]) for the potential V(x, y) from (a), Gaussian 
white noise £ x {t), and a dichotomous noise £ y (t) that flips at 
a rate 7 between ±A. Using dimensionless units, the param- 
eters values are r\ = 1, L = 1, B = 1, 6 = 0.55, 9 = 45°, 
ksT — 0.1, A = 10, 7 = 0.4. Dots with error bars: numerical 
simulations of (Q) . Solid line: Analytic approximation ([|)-(|^). 




FIG. 2. Hard-wall potential V(x,y) like in Fig. la but 
now periodically continued along the z-axis, resulting in a 
2-dim. array of obstacles ("sieve"). For symmetry reasons, 
the {/-component of the the Brownian motion (|l|) is completely 
independent of whether the grey shaded "corridor" is endowed 
with perfectly reflecting "side- walls" (Fig. la) or not. 
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FIG. 3. The current (El) versus the static force F in t\w, 
(g) for the potential V(x, y) from Fig. la or Fig. 2, indepen- 
dent Gaussian white noises £ x (t) and £th(t), and a periodic 
driving f(t) that jumps between ±A with time constant r ac . 
The parameters values are L — 4/xm, B — 3 /im, b = 1.2 fim, 
9 = 70°, T = 293° K, A = 0.2pN. The Brownian particle is 
assumed as spherical with radius r = 0.5 \xm and as subjected 
to Stokes friction r\ — iS-Kvr, where v is the viscosity of water. 
Dots with error bars: numerical simulations of (Q), with 
T ac = Is and Tac = 25 s in (a) and with r ac — 5 s in (b). Solid 
line in (a) with r ac = Is: Analytic approximation @-(P) 
with p(t) = S(t — T ac )- Solid line in (a) with r ac = 25 s: An- 
alytic approximation (|^), (Q)-(^|) with r = r ac . Solid line in 
(b): Analytic approximation by adding up (^) and (Q) and 
then evaluating with (^|), (|^), p(r) = 8{t — Tac), and 
Tac — 5s. Dashed line in (b): Analytic approximation for 
Tac = 5s based on jic|), ([ll[), as described in more detail in 
the main text. Note that all analytical result are based on 
the assumption of a point particle. The finite particle radius 
r has been approximately accounted for by replacing B by 
B-2r in (i). 
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